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, Abstract 

00 ■ We prove a regularity result for solutions of a purely integro-differential Bellman equation. 

This regularity is enough for the solutions to be understood in the classical sense. If we let 
the order of the equation approach two, we recover the theorem of Evans and Krylov about 

< j the resu '" i,y ol soluti °™ to — 8q " ation " 

^ . 1 Introduction 

Cq . 

In 1982, L. Evans and N. Krylov proved independently ([6] and [7]) the following celebrated interior 
regularity result for elliptic partial differential equations: If u is a bounded solution to F{D 2 u) 
\ in Bi, where F is uniformly elliptic and concave, then u € C 2 a (Bx/2) for some a > 0. In this 

■ paper we prove a nonlocal version of that theorem. We prove that solutions to concave integro- 
CTN , differential equations of order a have regularity C a+a for some a > 0. This is enough regularity 

to consider the solutions to be classical. 

The equations we study arise in stochastic control problems with jump processes (see for 
example [9], [8]). In [9] a C 2,a regularity of the solutions of Bellman equations for Levy processes 
is obtained, but the equation is required to have a uniformly elliptic second order part which 
is ultimately the source of the regularity. In pQ a purely integro-differential Bellman equation 
is studied. They only consider the case of the maximum of two linear operators. They obtain 
solutions in the fractional Sobolev space H a l 2 up to the boundary and H a in the interior of the 
. ! domain. As they point out in their paper, the solutions to these equations are expected to be 

■ more regular in the interior of the domain. 

In this paper we consider purely integro-differential equations and obtain an interior regularity 
result. Since we do not require our equations to have a second order part, our estimate comes 
only from the regularization effects of the integrals. 

The constants in our estimates do not blow up as a — > 2, so we can recover the usual Evans- 
Krylov theorem as a limit case. It is interesting to follow what the ideas of the proofs become as 
(j — ► 2. Interestingly, the ideas we present in this paper provide a different proof of the Evans- 
Krylov theorem for second order elliptic equations. 

We consider the equation 

Iu{x) := inf L a u{x) = inf / (u(x + y) + u(x — y) — 2u{x))K a {y) dy = (1-1) 

As in [3], we will choose each linear operator L a in some class C. Consequently, the operator / 
will be elliptic with respect to C in the sense described in [2]. 

We describe below the appropriate classes of linear operators that we will use in this paper. 
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We say that an operator L belongs to Cq if its corresponding kernel K satisfies the uniform 
ellipticity assumption. 

The ellipticity assumption (|1.2p is the essential assumption that leads to a local regularization. 
Our proofs, as usual, involve an improvement of oscillation of the solution to the equation (or an 
operator applied to it) in a decreasing sequence of balls around a point in the domain. Since the 
equations are nonlocal, every argument in our proofs will have to take into account the influence 
of the values of the solution at points outside those balls. We will often need to say that the part 
of the integral in (jl.ip outside a neighborhood of the origin is a smooth enough function. That is 
why we define the following classes of smooth kernels. 

We say that L g L\ if, in addition to Q1.2p . the kernel K is C 1 away from the origin and 
satisfies 

VK a (y) < (1.3) 
Finally, we say L £ £2 if the kernel is C 2 away from the origin and satisfies 

D 2 K a iy) < , , ^ (1.4) 

We consider the corresponding maximal operators 

M+u{x) 
M^u(x) 



sup Lu(x) 
LeCo 

sup Lu{x) 
LeCi 

sup Lu(x) 

LGC 2 



A(5u{x,y))+ - \(Su(x,y)Y 



\n+a 



dy 



Recall that we write 5u(x, y) = {u{x + y) + u(x — y) — 2u(x)) as in .3;. Naturally we have the 
inequalities M^u > M^u > M^u. The minimal operators M~ are defined likewise. 

We do not know a closed form for or M + 2. Since £2 C £1 C £0, we have the relations 
M%u < Mfu < M^u and M^u > Mfu > M$u. 

Our main result states that under the hypothesis that all operators L a belong to £2, the 
solutions are classical in the sense that there is enough regularity so that all integrals are well 
defined and Holder continuous. 

Theorem 1.1. Assume every L a in (|1.1|) belongs to the class £2. If u is a bounded function in 
R" such that Iu = in B 1} then u £ C a+a (B 1/2 ). Moreover 

\\ u \\c°+«(B 1/2 ) < C \\ u \\l°°(R«) ( L5 ) 

The interested reader may verify, in following the arguments, that global boundedness of u 
may be substituted by an appropriate moderated growth at infinity (see remark at the end of this 
paper). 

For values of a less or equal to 1, this theorem does not provide any improvement with respect 
to the C 1,Q estimates in [2]. Thus, for the purpose of proving Theorem ll.il we will assume a to 
be strictly larger than 1 in this paper. The result becomes most interesting when a is close to 2 
and a + a > 2. 
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Note that the result of the theorem remains true if I is convex instead of concave (a sup 
of linear operators instead of an inf). Indeed, we can transform one situation in the other by 
considering the equation —I(—u) = 0. 

In previous papers ([3J and [5]) we started developing the regularity theory for nonlocal equa- 
tions. In [3J we obtained a nonlocal version of Krylov-Safonov theory with estimates that do not 
blow up as a — > 2. This allowed us to obtain C 1,a estimates for general fully nonlinear integro- 
diffcrential equations that are translation invariant. In [2], we extended those results to variable 
coefficient equations using perturbative methods. In this paper, we use the results in our previous 
two papers extensively. 

2 A regularization procedure 

In this section we show a simple technique to approximate uniformly the solutions to the integro- 
differential equation Ql.ip by C 2,a functions that solve an approximate equation with the same 
structure. This procedure works exclusively for integro-differential equations and cannot be done 
using only second order equations. It makes it unnecessary to use sup- or inf- convolutions and 
simplifies the technicalities of several proofs. Essentially the idea is that if we prove an estimate 
assuming the solutions u is C 2a (but the estimate does not depend on the C 2,a norm), then we 
can pass to the limit using this approximation technique to extend the estimate to all viscosity 
solutions. In this respect, the technicalities in the integro-differential setting simplify very much 
compared to the second order counterpart. 

Lemma 2.1. Let u be a continuous junction in W 1 solving (jl.ip with every L a belongs to the class 
£2 (resp. C\ or Co). There is a sequence of regularized equations in the same class 



so that the solutions u e are C 2a in the interior of B\ for every e > and lim £ ^o u £ = u uniformly 
in B\. 

Proof. Let r\ be a smooth function such that 



I e u £ = inf L E a u £ = 



in B\ 



a 



u e — u 



in W 1 \ B x 



< 77 < 1 in K™, 
77 = in R" \ B 
7] = 1 in B\ , 



and let r\ e (x) = r/(x/e). 

Let us consider the following regularized kernels 



2-0- 



+ (1 - r, e {x))K a {y) 



\v\ 



Correspondingly we define 




I s v = inf L e a 



v 



Note that if L a G C{ then also L £ a £ Li for i = 0, 1, 2. 



3 



Let u e be the solution of the following Dirichlet problem: 



I £ u £ = ML £ a u £ = 



in B x 



a 



u £ = u 



in E n \ B x 



The solution u £ to this problem is C 2,a by theorem 6.6 in [2J. 

It is clear that if v € C 2 (x) and \v(y) — v{x) — (y — x) ■ Vv(x)\ < M\y — x\ 2 in B-y, then 



\Fv(x) -Iv(x)\ < CMs 



so ||/ e — 1\\ — > as e — > (recall cr < 2), where the norm ||/ e — I\\ is computed in the sense of 
definition 2.2 in [2J. Then, by lemma 4.9 in [2J, u e converges to u uniformly in B\ as e — > 0. □ 

Remark 2.2. The concavity of / is not used in Lemma 12.11 The exact same idea works for 
equations of the type 



3 Average of subsolutions is a subsolution 

The main ingredient in the Evans-Krylov theorem is the fact that concavity of the equation makes 
second order incremental quotients subsolutions of the linearized equation. In order to prove that, 
one first observes that an average of solutions to a concave equation is a subsolution to the same 
equation. 

In this section we prove that also in the non local case the average of subsolutions to a concave 
equation is a subsolution of the same equation. This is obvious from the equation if the solutions 
are classical. For viscosity solutions we can prove it quickly using the approximation technique of 
section [5] 

Proposition 3.1. Let u and v be subsolutions of Iu = and Iv — in a domain fl, u, v continuous 
in W\ then I(u + v)/2 > in Q. 

Proof. The proposition is obvious if it, v £ C 2 by the concavity of /. So we used the regularization 
procedure described in section [2] 

Let I £ u £ — and I £ v £ = be the approximate equations of Lemma 1 2. II The functions u £ and 
v £ are C 2 so I £ (u £ + v £ )/2 > in f2. Since u £ — > u and v £ — > v uniformly in fl and I £ — » /, then 
I(u + v)/2 > in Q by Lemma 4.9 in [2j. □ 

The same idea shows that any average of solutions is a subsolution. In particular we have the 
following proposition. 

Proposition 3.2. Let u be a solution of Iu = in B\ and rj be a mollifier, i.e. 

1. n>0 

2. fr) = l. 

3. supp r\ C Bs ■ 

we have I(n * u) > in Bis. 
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4 The linear theory of integro-differential equations 



In this section we present some regularity theorems for linear integro-differential equations with 
constant coefficients. Naturally in this simple case, we can easily obtain more powerful results 
than for the nonlinear case. The results we present in this section are just the ones that we will 
need in the rest of this paper. 



Theorem 4.1. Let L be an integro-differential operator in the class C\ with a > <7q > 1. Suppose 
that u is an integrable function in the weighted space L 1 (M n , ppj^p ) that solves the equation Lu = 
in B\, then u 6 C 2 ' a (Bi/ 2 ) and we have the estimates 

IMI c2 , Q(Si/2) <c fe |H| L1(R „ i j 



The value of the constant C and a depends on n, A, A and <7o but not on a. 

Proof. We will prove the apriori estimate. The regularity estimate for a weak or viscosity solution 
follows by mollifying the solution or using the regularization procedure of the previous section. 
First we apply theorem 2.8 in [2] to obtain that u £ C 1 ' a (B 3 / 4: ) and obtain the estimate 

\\u\\cHB 3/s ) <C\\u\\ Ll{R ^_^ ) . 

The idea is to apply the same C 1,a estimate to every directional derivative u e . Since we do 
not have an L°° estimate of u e outside of -B3/4 we have to use our usual integration by parts trick. 
We know that 

u e (y)K(x + y) dy = 



Let 77 be a smooth cutoff function such that 

< r] < 1 
r] = 
rj = 1 



m K 

outside B- 
in B 5/s . 



3/4, 



We compute 



u e (y)v(y) K (x + y) dy 



u e (y)(v(y) - l)K(x + y) dy 



Ue(y)(Ve(y)K(x + y) + (r)(y) - l)K e (x + y)) dy 



<C\\u\\ 



LH 



Thus we can apply Theorem 6.1 in [2] to conclude that u e £ C 1:Q for every direction e and 
thus it £ C 2,a (Note that we are using a > 1 here). □ 

In order to have better interior regularity estimates than C 2 ' Q , we would need to impose more 
regularity to the kernel K in L than C 1 away from the origin. 

Next theorem says that in L 2 all linear operators have a comparable norm. 

Theorem 4.2. Let Lq and L\ be two linear integro-differential operators in the class Co- Suppose 
that L a u £ L 2 (R n ) then L x u £ L 2 (R n ). 
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Proof. Since we are dealing with L 2 norms, and translation invariant linear operators, we will use 
the Fourier transform to prove this theorem. 

Given a function u and y £ R™, we have Su(x, y) = (u(. + y) + u(. —y) — 2uJ~~ = [e iy '^ + e~ ly '^ — 
2)"(£) = 2(cos(y ■ — !)"(£)■ We use this identity to compute the symbol s(£) of an operator —L 
as a pseudo differential operator. 

-MO = [J 6u(x,y)K(y) dy^) 

= (J 2(1 - cos(y • 0)K(y) dy) u(0 =: 

Note that for every £ function (1 — cos(y • £)) is C 2 an( i bounded, so the integral in the right 
hand side is well defined. Let us estimate it from above and below. 
For any R > 0, 

s(0 = [ 2(1 - co S (y ■ 0)K(y) dy < [ 2\y ■ £| 2 (2 - a)-^ dy + f 2(2 - a)-^ dy 

jr" Jb r w\ J» n \B B \y\ 

< C\£,\ 2 R 2 -° + C^— ^R- a 
a 

so we obtain s(£) < C\£\ a by choosing R = \£\~ x ■ 

On the other hand note that (1 — cos(y ■ £)) is nonnegative and so is K (y). So the integrand 
is nonnegative and we have 

s(£) = I 2(1 - cos(y • 0)K{ V ) dy > f -\y ■ £| 2 (2 - a)^^ dy 

So the symbol s(£) is comparable to |£| -<T for any operator L in Cq. Thus by classical Fourier 
analysis we have that LiLq 1 has a bounded symbol and maps L? functions into L 2 . □ 

The following theorem is a direct combination of Theorems 14.11 and 14.21 

Theorem 4.3. Let L be an integro- differential operator in the class C\ with a > ctq > 1. Suppose 
that u is a function in the weighted space L 1 (M n , ppj^p ) that solves the equation Lu — f in Bi 
for some f £ L 2 . Let L\ be an operator in Cq, then L\u £ L 2 (i?!/ 2 ) and 

IM| i2(Bi/2) < c (||u|| L i (R „ i _i_ )+ || / | L2(bi J 

for some constant C depending on n, A, A and a®. 
Proof. Consider the function v that solves 

Lv = f XBl in R". 

From Theorem 14.21 we get that L\v £ L 2 (R"). Theorem 14.21 can also be applied to L^ — 
(-A) CT / 2 , thus v is in the homogeneous fractional Sobolev space IP 7 / 2 . By Sobolev embedding 
v £ LP(W l ) where p = An/(2n - a). In particular v £ L 1 (W l , j+j^p) 

Now we apply theorem 14. II to u — v £ L 1 (M. n , pp^p ) to conclude the proof. □ 

Remark 4.4. The cllipticity constants A and A of Lq and L\ do not need to coincide. Indeed if 
each Li is elliptic with constants A, and Aj, then they would both be elliptic with respect to the 
constants min(Ao,Ai) and max(Ao, Ai). 
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5 Subsolutions in L 1 are bounded above 

The theorem below is a weak version of the mean value theorem. Its proof uses the same ideas as 
the proof of the Harnack inequality in [3]. We include it here for completeness. 

Theorem 5.1. Let u be a function such that u is continuous in B\, assume that 

- H f + dy<C 

i + \y\ n+a ~ 

M+u > -C in Bi 

then 

u(x) < CCq in B1/2 
for every x G B1/2, where C is a universal constant. 

Proof. Dividing u by Co , we can assume without loss of generality that Cq = 1 . 
Let us consider the minimum value of t such that 

u(x) < h t (x) := t(l - \x\)~ n for every x G B x . 

There must be an xq G B\ such that u(xq) = ht{xo), otherwise we could make t smaller. Let 
d = (1 — \xq I) be the distance from xq to dB\. 

For r — d/2, we want to estimate the portion of the ball B r (xo) covered by {u < u(xq)/2} and 
by {u > u(xq)/2}. We will show that t cannot be too large. In this way we obtain the result of 
the theorem, since the upper bound t < C implies that u(x) < C(l — 

Let us first consider A := {u > u(xq)/2}. By assumption, we have u G thus 

\AnBx\ < c 



u(x ) 

< cr x d n 

Whereas \B r \ = Cd'\ so if t is large, A can cover only a small portion of B r (xo) at most. 

\{u > u{x Q )/2} n B r (x )\ < Cr 1 \B r \ (5.1) 

In order to get a contradiction, we will show that \{u < u(xo)/2} n B r (xo)\ < (1 — a)B r for a 
positive constant a independent of t. 

We estimate \{u < u(xo)/2} n Be r (xo)\ for > small. For every x G Bg r (xo) we have 
u(x) < h t (x) <(d- 8d/2)- n < u(x )(l - 6/2)~ n , with (1 - 0/2)-" close to one. 

Let us consider 

v(x) = (1 - 6/2)- n u{x Q ) - u(x) 

so that v > in Bg r , and also M~v < 1 since M + u > — 1. We would want to apply Theorem 10.4 
in [3] (the L e estimate) to v. The only problem is that v is not positive in the whole domain but 
only on Bg r . In order to apply such theorem we have to consider w = v + instead, and estimate 
the change in the right hand side due to the truncation error. 

We want to find an upper bound for M~w = M~v + instead of M~v. We know that 

M-v(x) = (2-a) [ ^(x, g )+ - Mt,(x,y)- 

jR n \y\ 
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Therefore, if x S Bg r / 2 (xo)i 



X5w(x, y) + — ASw(x, y) 



M-w = (2 -a) v tttz V ' dy (5.2) 



v(x + y) 



< 1 + ( 2 - °) / -A l An+ l' dy (5.3) 



i n{v(x+y)<0} \U 



< 1 + (2 - a) [ A (u(x + y)-(l e/2)-Mxo)r dx (g>4) 

jR n \B er/ 2 \y\ 

<l + C(2-a)(er)- n - a f A—^p^— dx <c(dr)- n - a (5.5) 



Now we can apply Theorem 10.4 from [3] to w in Bg r / 2 {xo)- Recall w(xq) = ((1 — 0/2)~ 
1)u(xq), we have 



u < 



Y^j n % (xo) = |{«J > u(a;o)((l - 0/2)- - 1/2)} n B 9r/i (x )\ 

< C(6r) n (((l-6/2)- n - l)u(xo) + C(6r)- n -° (r0) CT ) e (u(x )((l - 0/2)-" - ~) 

< C(0r)" (((1-0/2)-" - 1) £ + 0-" £ t- £ ) 

Now let us choose > so that the first term is small: 

C(0rr((l-0/2)-« -If <\\B 6r/2 \ . 

Notice that the choice of is independent of t. For this fixed value of we observe that if t is 
large enough, we will also have 

C(0r)"0-" £ i- £ <\\B 9r/2 \ 

and therefore 

|{it < w(x )/2} n Be/r4(^o)| < 2 | B er/4(a;o)| 
which implies that for i large 

|{u > «(a;o)/2}n4/4(j;o)| > c\B r \ . 
But this contradicts (|5.1|) . Therefore i cannot be large. Rescaling back, we obtain 

u{x) < CC , 

for any x in B\j 2 . O 

6 Each L a w is bounded 

The idea of this section is to show that averages of second order incremental quotients are subso- 
lutions (of the maximal operator M + ). In particular, each v a := L a u is a subsolution to 

Va > £ Bi 

M + v a >-CeB l 

Then we would estimate the integral of v a in £>i/2 and use Theorem 15.11 to prove that L a u is 
bounded. 
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Lemma 6.1. Assume u e C 2 . Then for every kernel K corresponding to an operator L £ Cq and 
every bump function b such that 

< b(x) < 1 in R n , 
b{x) = b(-x) in R™, 
b(x) =0 in M" \B 1 / 2 , 



have 



M+ Qf Su(x,y)K(y)b(y) dy^) > m B 1/2 



Proof. Let 0fc be the L 1 function 4>k{y) = XR n \B 1/k (y)K(y)b(y). Since u 6 C 2 , we can approximate 
the value of the integral uniformly by 

Su(x, y)K(y)b(y) Ay = lim / Su(x,y)(f>k(y) dy = lim u * 0* - ||^fe|| L i u. 

Applying Proposition ^. 11 we have 

'("rar)*° 

On the other hand, we know that Iu = 0, then by Lemma 5.8 in [3] and the fact that M + is 
homogeneous, 

•/. 



M+(u*0 fe - Uk\\ L iu) = ||0 fc || Ll M+ (u* 



-u > 



Since we have M 2 (u * (j)f. — ||</>fc|| L i u) > for every fc > 0. Then the result follows by Lemma 
4.3 in [5] by taking limit as k — ► oo. □ 

Lemma 6.2. Assume u 6 C 2 . TTien i/iere is a universal constant C such that for every operator 
LeC 2 

M+ (Lu) > -C in B 1/2 

Proof. As in the proof of Lemma l6.1i we let <f>k be the L 1 function 4>k{y) — XR n \B 1/k (y)K(y) and 
we approximate the value of the integral uniformly by 



Lu(x) = lim / Su(x,y)<f> k (y) dy = lim (u * 4> k - u \\4>k\\ L , ) 

k — >oo J^n k — >oo 

Let & be a bump function such that 

< b(x) < 1 in K", 
b( x ) = b(-x) in R", 
b{x) = in W l \B 1/2 , 



As in the proof of Lemma I67T 



Therefore (u * (0 fe &) - ||0fc&|| L i u) > 
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On the other hand, we estimate I(u * ((j) k (l — &))) in Bi/2- 

I(u * (cf> k (l - b))) = inf / L a {u * (0 fc (l - 6))) chy 
= inf / u * L a ((f) k (l - b)) dy 

a JR™ 

Since by fL~2]l . </> fc (l - 6) e i 1 and by |[L~3]) . D 2 4> k (l - b) G L 1 uniformly in k, then L{cj) k (l - b)) 
is bounded in L 1 uniformly for all L G Cq. Therefore, 

\I(u*(cf> k (l-b)))\<C\\u\\ Lco . 

Using Lemma 5.8 in [3] and the homogeneity of M~, 

M 2 - (u * (0 fc (l - &)) - ||0 fe (l - u) > -C \\u\\ Lao 

Therefore 

M+ (u * <f> k - \\<p k \\ Ll u) >M+(u* {cj> k b) ~ U k b\\ Ll u) + M- (u * (<fc(l - 6)) - ||0 fc (l - 6)|| LI u) 
> -C\\u\\ LBB 

We finish the proof of the lemma by taking k — > oo using Lemma 4.3 in [2]. □ 

Lemma 6.3. Let u be a solution of Assume L a G £2 /or every a. Then for every a, 

L a u < C \\u\\ LOO in -B1/8 for some universal constant C. 

Proof. By Lemma 12.11 we can assume the function u is C 2 . Indeed, for every e > we can 
approximate u with a C 2 function u s that satisfies the same kind of equation. If we can prove the 
estimate for u e with a universal constant C that does not depend on e, then we would prove it for 
u by passing to the limit as e — >■ 0. So we assume that u G C 2 and thus all the integrals are well 
defined. 

From Lemma 16.21 we know that for each L a , 

M+(L a u) > -C\\u\\ Lac mB 1/2 . 

We would want to apply Theorem 15.11 to L a u. For that we still need an estimate at least in 
L 1 ((l + |y|)~" -cr ). We can easily obtain an estimate in L 1 (B 1 / 2 ) using the fact that L a u > in 
Bi because of the equation (jl.ip . 

Let b be a smooth cutoff function such that < 6 < 1 in W 1 , b = 1 in -B1/2 and 6 = outside 
B\. We multiply L a u by b and integrate by parts. 

I L a u{x) b(x) dx < L a b(x) u(x) dx < C\\u\\ Lao 

for some universal constant C. Since L a u > in B\, then it is in L 1 (B 1 / 2 ). 

We would still need some control on the values of L a u away from B 1 / 2 in order to apply 
theorem 15.11 We do not want to assume any regularity for u outside B\ , so our only choice is to 
cut off again. 
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Let c(x) := b(2x) and w(x) = c(x) L a u(x) We will estimate w(x) for x S B1/4. For that, 
let us consider any operator L G £2 and estimate 



Lw{x) = j 8w(x,y)K(y) dy 

= I 5L a u(x,y)K(y) dy S(L a u(l - c))(x, y)K{y) dy 

JR™ JR™ 

> -C ||«|| ioo - 2 / L a u(x + y)(l- c(x + y))K{y) dy 

> -C \\u\\ LOO -2 J u(x + y)L a ((1 - c(x + .))K) dy 
>-C\\u\\ L - 

For the last inequality it was used that L a ((1 — c(x + -))K) is in L 1 uniformly for x S B1/4. 
This follows from the estimates (|1.2[) and (| 1 .4|) of the kernel K. 

Now we can apply Theorem 15.11 to w to obtain that w < C ||u|j Loo in Bi/ S . but w = L a u in 
i?!/4, so we finish the proof. □ 

7 Extremal operators are bounded 

In this section we prove Theorem 17.41 We will achieve this result by showing that Lu(x) is 
bounded uniformly for all L 6 Co- The ideas are very similar to section [6] but now we apply 
them to operators that are not a priori bounded below, so we use the I? estimates from section 
Q] instead. 

Lemma 7.1. Assume u 6 C 2 . T/ien /or every kernel K corresponding to an operator L £ Co and 
every bump function b such that 

< b(x) < 1 m R", 
b(x) = b(-x) in R n , 
b(x) = m R n \B 1/2 , 

then 

M+[b(x)[ 5u(x,y)K(y)dy)>-C\\u\\ Lao in B 1/2 

\ JB 1/2 J 

Proof. Let us call 

L t u(x)= / Su(x,y)K(y) dy. 



J B1/2 

r+rrt. 



By Lemma RTT1 we have Mj(L'u) > 0. Let L be any operator in £2, so we estimate 
L(bL t u)(x)= f 8{L t u){x,y)K{y)dy- f S((l - b)L t u)(x 1 y)K(y) dy 

JR" JR™ 

> -2 j£ (l - 6(x + y))^^ + y)A'(j/) dy 

> -2 y ^ u(a; + y)L t ((l - 6(1 + dy 
>-C\\u\\ LOO 

where we used that — 6(x + -))K) is bounded in L 1 uniformly in x. This is due to the fact 

that L> 2 ((1 - b(x + -))K) e L 1 (R n ) because of {HI- □ 
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Lemma 7.2. Let u be a solution of with all operators L a in £ 2 . There is a constant C such 

that for every operator L in Cq we have 

\Lu(x)\ <C\\u\\ LOO infli/a 

Proof. As in the proof of Lemma l6.31 we can and will assume that u e C 2 . We will write the proof 
assuming that H^H^^n) = 1- Moreover, we will prove the estimate in Bi/64 instead of Bi/ 2 . The 
general estimate follows directly by scaling and a standard covering argument. 

Let L a be one of the operators used in the infimum in (jl.ip . We know from Lemma 16.31 that 
L a is bounded in B 1 / 2 by a constant C. In particular ||-kaw|| i2 ( Sl/2 ) < C. From Theorem 14.31 we 

have an L 2 estimate for every L £ Cq, 

\\Lu\\ L2{Bi/i) < a 

We split the integral of Lu into two domains 

Lu(x) = / Su(x,y)K(y) dy + / Su(x,y)K(y) dy. 

The second integral is clearly bounded since K(y) is a function in i 1 (R™ \B 1 / 2 )- Thus we still 
have an estimate in L 2 for the first term. 

11 < C 

L 2 (B 1/4 ) 



Su{x,y)K(y) dy 

2 

On the other hand, from Lemma |6.1[ 



M+ / Su(x,y)K(y)dy) >0 

\JB 1/2 ) 

For a bump function c(x) such that 

SUppC = 

c = 1 in B-i/g 

we define 



w(x) :— c(x) I Su(x,y)K(y) dy 

We know that w 6 L 2 (R") and w — outside B 1 / i . In particular w is bounded in the weighted 
L 1 space: L x ((l + |y| n+ff )~ 1 ) needed for TheoremO 

We estimate M+u> > -C in / 16 now in the same way as in the proof of Lemma 16.31 (using 
assumption (|1.4[1 ). Thus we can apply theorem 15.11 to w to obtain that w < C in -B1/32, which 
naturally implies L k u < C in Bi/ 32 - Passing to the limit we get that Lu < C for any L £ Cq. 

We got the desired bound from above only. In order to get the corresponding bound from 
below we must use the equation. 

Recall that L a u is bounded by Lemma 16.31 and the formulas of L a and L are given by 



L a u(x) = j 5u(x,y)K a (y) dy 
Lu(x) = / 8u(x,y)K(y) dy 
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where both K and K a are bounded below by (2 — s)\/\y\ n+a and above by (2 — s)A/\y\ n+<T . 
Consider the kernel 

2 1 
K d = -K a - —K 
A A 

and the corresponding linear operator L d . The kernel K d satisfies the ellipticity conditions (2 — 
v)/\y\ n+a <K d <(2~ cr)(2A/A - X/A)/\y\ n+a , so L d is in the class C with ellipticity constants 
1 and (2 A/ A — A/ A). The same proof as above tells us that L d u < C in Bi/ 32 . But then since L a 
is bounded, then we obtain a bound below for L in B\ / 32 . 

Lu = 2^L a - AL d > -C 
A 



Thus, we have both bounds and we obtain | \\Lu\\ \ < C in Bi/ 32 . 
Corollary 7.3. M^u and M^u are bounded in B 1 / 2 . 



□ 



Proof. Since M^u = sup ig £ Lu and for every L in Co we have \Lu\ < C \\u\\ Lxi with C indepen- 
dent of the choice of L in Co, then also |Mq~u| < C \\u\\ LO o with the same constant C. □ 

Now we are going to prove that all integrals in are absolutely convergent. This already 
implies that the solution is classical in some way. 

Theorem 7.4. Assume every L a in belongs to the class C 2 . If u is a bounded function in 

1™ such that lu = in B\ in the viscosity sense, then we have the following estimate 

[ \5u(x,y)\ |f dy < 0\\u\\ Loa(&n) in B 1/2 

jR n \y\ 

Proof. Applying Lemma I7T21 to L = — (— A) CT / 2 we get 
\-(-A)°/'u(x)\ = 
On the other hand, applying Corollary [731 with any pair A < A 



Su(x,y)j—^ dy 



<C\\u\\ LOO(Rn) inB 1/2 . 



\M+u(x)\ = 
Subtracting, we obtain 



(2 - s) 

{A6u{x,y)+ - X6u(x,y) )±—^-dy 



<C\ 



in B 



1/2- 



M+m(x) + A(-A) <T/2 u(x) < C ||«|| io 
(A -A) / 5u (x,y) + ^-^-dy<C\\u\\ LO 



On the other hand, by subtracting A(~A) a / 2 u(x) — Mqu(x) we obtain the bound 

(A - A) J 5u(x, y) _ fcj$ dy = A{-Ay/ 2 u{x) ~ M+u(z) < C \\u\\ L ^ n) 
Combining the two estimates above, we finish the proof. 



□ 
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8 Outline of the strategy: the second order case. 



Theorem 17.41 provides an estimate slightly stronger than u E C a . In the case cr — ^ 2 it becomes an 
estimate of the C ' norm of u. Comparing with the proof of Evans- Krylov theorem (as in [6], [7] 
or [1]), the underlying strategy of the proof up to this point is essentially the same but adapted 
to the integro-differential setting using the ideas in our previous papers [3] and [3J. 

The next step in the proof of our main result is to pass from this C a estimate to a C a+a esti- 
mate. In the second order case it corresponds to the apriori estimate ||w|| C 2.a,- Bi/2 ) < C 
The presently known proofs of this apriori estimate seem difficult to adapt to the nonlocal setting. 
Thus we present a different strategy for the proof. The key tools that the proof is based on are 
similar, but in our approach they are organized differently and arguably more directly. We plan to 
publish a short note [5] focusing only on this new proof for concave second order elliptic equations. 

In order to better understand our proof in the next section, we first sketch its adaptation to 
the second order case. 

We consider a C 2 solution of a fully nonlinear equation F(D 2 u) = with F concave and 
uniformly elliptic. We transform this into an integral equation by pointing out that a linear 
equation can be written as an integral on the unit sphere S\. 



a,ijdiju{x) = I d aa u(x) w(a) da 

with the weight w(a) = l/(det{aij}a l: >aiaj) where {a 1 - 7 } = {ctij} -1 . If the coefficients a^- are 
uniformly elliptic, then w(a) will be bounded away from zero. 

We recall that F being concave implies that pure second derivatives are all subsolutions of the 
linearized operator. In particular, for any fixed set A C Si, so is 

va = d aa u{x) dcr. 
J A 

We also recall that for a (non necessarily concave) fully nonlinear equation 

F(D 2 u) = supinf L ab u = 0, 
b a 

a solution u satisfies (just because it is an inf sup) that for any two points x and y in the domain, 
there exists an operator L a b for which 

L ab u(x) - L ab u(y) > 0. 

In our approach, this means that there is a weight w(a), bounded below and above depending on 
the ellipticity constants, such that 



(d aa u(x) ~ d aa u(y)) w(a) da > 0. 



Si 

In particular, since there is another weight which gives the same inequality exchanging x and y, 
we must have that the following quantities are comparable: 

{d aa u{x) - d aa u{y)) + da k I (d aa u(x) - d aa u{y)y da « / \d aa u(x) - d aa u(y)\ da 

Si ^ Si ^ Si 



At this point we define 



(8.1) 



h(x, a) = d a(T u(x) - d aa u(0) 
wa(x) = / h(x,a) da 

J A 
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for any set A C Si. We will use only the properties above to show that 



/ \h(x,a)\ da < \x\ c 



The C 2,a estimate for u follows easily from this estimate. 

By (|8.1[) . we only need to prove that wa(x) < C\x\ a for every set A, since 



\h(x,a)\d~ / h(x, ay da = sup wa(x) 

JS! A 

In fact, by renormalization we only need to prove the following lemma. 

Lemma 8.1. Assume that for x in B±, for any set A C Si, wa(x) < 1. Then there is a universal 
constant 8 > such that 

w A (x) < 1 - 6 

for any A C Si and x £ -Bi/2- 

Sketch of the proof. Suppose that there exists an x € -B1/2 where wa > 1 — for some set i C Si. 
We will arrive to a contradiction if 9 is too small. 

Since wa < 1 in Bi and wa is a subsolution of the linearized equation, we can apply Theorem 
4.8(1) in 3] (the L e estimate) to 1 — wa (this will correspond to Theorem 10.4 in [3] in the nonlocal 
case). It follows that we can make 

n = {w A {x) >i-te} 

cover almost all -B1/4 if we choose t large (but independently of 8). 
We will now obtain a contradiction by looking at wa<= in -B1/4. 
For every x in f2, the choice of the set A is almost maximal in the sense that 

1 - tO < w A {x) < / h(x,a) + da < 1 

On the other hand, since 

/ h(x, a) da = / h(x, a) + da — I h(x, a)~ da = wa(x) + wa<= (x) 

J Si •'Si ^oi 

then also 

< wa c + / h(x, a)^ da < t9 
From (|8.ip . we know that the integrals of h + and h~ are comparable. Thus in Q we have 

wa- <te-c 

for a constant C depending on A and A. If we choose small, that means that wa- will be strictly 
negative in most of -B1/4. But then applying Theorem 4.8(2) in [1] (which corresponds to Theorem 
15.11 in the nonlocal case) we obtain that wa-(0) < —c for some universal constant c. This is a 
contradiction since clearly ^^^(O) = 0. □ 

In the integro-differential case, the proof will be slightly lengthier in part because we have to 
keep track of the truncation error we make every time we localize an integral. We cover the proof 
in detail in the next section. 
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9 Further regularity. 

This section is devoted to fill the gap between Theorem 17.41 and Theorem 11.11 
From Theorem 17.41 we know that 

2 - (7 



V 

Our objective is to show that 



' \5u{x iy )\-—&y<C in B 1/4 (9.1) 

Bi/2 \y\ 



\Su(x,y)-Su(0,y)\^—^-dy<C\x\ a 



b 1/2 \y\ 

for some constant C and a > and for every x £ -B1/4. This estimate implies the Holder continuity 
of the fractional laplacian (— A) CT / 2 from which the C a+a regularity of u follows. 
We will consider all kernels K of the form 

where Xa(u) is the characteristic function of an arbitrary set A. 

Let b be a bump function as in Lemma l7.ll For each set A, we write 



w A (x) = b(x) / (Su(x, y) - Su(0, y))K A (y) dy 
- '-B1/2 

We know that wa is uniformly bounded from Lemma 17.21 From Lemma 17. 11 we have 

M^wa > -C in B 1/4 uniformly in A (9.2) 
We define the following quantities 

P(x) := supwa(i) = b(x) [ {5u{x 7 y) - Su(0, y)) + y °] dy 
a Jb 1/2 \y\ n+a 

N(x) := snp-w A (x) = b(x) [ (Su(x, y) - Su(0, y))~ 2 ° dy 

a J Bl/2 \y\ n+a 

Note that P(x) is realized by the set A = {x : 5u(x, y) > Su(0, y)} and N(x) is realized by the 
complement of that set. 

Lemma 9.1. There is a constant C such that for x £ B 1 / 4 , 

yJV(a?) -C\x\< P(x) < ^N(x) +C\x\. 
A A 

Proof. For some x £ -B1/4, let u x (z) := u(x+z). Since u solves the equation (|1.1[) in a neighborhood 
of x, then both u and u x solve (jl.ip in a neig hborhood of 0. Thus M+(u x - tt)(0) > and 
M2 (u x -u)(0) < 0. 

For every kernel K in the family £2 we have 



L{u x -u)(0)= I (Su(x,y)-Su(0,y))K(y)dy 
(6u(x, y) - Su(Q, y))K(y) dy 

(6u(x, y) - 6u(0, y))K(y) dy 



1/2 
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Let us analyze the second term in the right hand side. 

(Su(x, y) - Su(0, y))K(y) dy = Su(0, y) (K(y - i)xbj /2 (y - x) + K(y)xBf /2 (y - x)j dy 



< 



\ SU (0,V)\ |..|n+«r+l M ^ + 8 ML" / \ v \n+a d y 

"\Si/ 2+ ixi ifl Jb 1/2+m \b 1/2 \y\ 



< C\x\ 

Therefore, for every kernel K in the family £2, we have 

(Su(x, y) - 5u(0, y))K{y) dy > f (Su(x, y) - 5u(0, y))K{y) dy + C\x\ 

JB 1/2 

Taking the supremum we obtain 

< M$(u x - u) < sup f (5u(x,y)-5u(0,y))K(y)dy + C\x\ 

K JB 1/2 

In particular, if we take the suppremum over all kernels K in Cq (a larger family), we still have 



I Ml 



sup / (Su(x, y) - 6u(0, y))K(y) dy > -C\x\ 

<K<K^z2^JB 1/2 



% +" — — 



which is the same as AP(x) — XN(x) > —C\x\. 

The same computation with M^" (u x — u) (0) < provides the other inequality. □ 

It is important to notice the following relation, 

2 - a 

\Su(x,y) - 6u(0,y)\ dy = supw A - Mw A 



v \n+o 

1/2 i y 



P(x) + N{x) 



The strategy for proving our regularity result will be to prove that sup^g^ P(x) < Cr a . It is 
enough to prove it for \x\ small enough, therefore we can consider a rescaled situation by taking 
wa(x) = ^jWA(rx) where C is the constant from (|9.ip and r is small enough so that our estimates 
become 

for every K: \w K \<l in R" (9.3) 
for every K: M+w K > -ei in B 1 (9.4) 

^N(x) - eiM 1 " 61 < P(x) < Y N(x) + eilxl 1 " 61 (9.5) 
A A 

for ei arbitrarily small. 

Lemma 9.2. Assume a € (T2). Let P(x) be the function defined above. There is a constant C 
and a > such that 

P(x) < Cr a 

Proof. As mentioned above, after an appropriate scaling, we can assume that (|9.3p . (|9.4|) and 
(|9.5p hold with e\ arbitrarily small. On the other hand, given the construction in Lemma [2~T1 we 
can assume u is C 2 and thus wk , P and A?" are continuous. We will obtain the apriori estimates 
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independently of the modulus of continuity of them, so the estimate holds when passing to the 
limit. 

We will prove that there is r > and 8 > such that 

sup |P| < (1 - 8) k = r ak where a = ^ - ^ (9.6) 
B rk logr 

This is clear for k — 0. Let us prove it is true for all values of k by induction. So let us assume 
it is true up to some value k. 

Since (|9.6p holds up to some value k, we have that, 

\w A {x)\ < (1 - 6»)- 1 |x| Q for |x| > r fc 

Consider the following rescaled functions 

w A (x) = (l-9)- k w A (r k x) 
P(x) = (1 - 6)~ k P(r k x) = supw A (x) 

A 

JV(a;) = (1 - ey k N(r k x) = sup-w A (x) 

A 

The function P satisfies the relations 

P(x) < 1 in B x 

P{x) < (1 - 0) -1 |a:| a outside Pi 

Moreover, from (|9.5p . 

A7V(x) - £i < P(x) < A7V(a;) + e x (9.7) 

We want to show that if 8 and r are chosen small enough we will have P < (1 — 6) in B r . The 
proof is by contradiction. We will arrive to a contradiction if 9 and r are small enough. 

Let xo be the point where the maximum of P is achieved in B r for some r 6 (0, 1/2). We 
assume P(xo) > 1 — to get a contradiction. Let A be the set such that P(xq) — wa{xo) >l — 0. 

Let va = (1 — We know that inf j B r > 6*. Moreover, 

MJv^ < MJ(1 - w A ) -MJ(1 - w A )~ 

< -M%wa - Mj (1 - 
<C inP 1/2 

since M+w A > si and (1 - w A y < ((1 - e)- 1 ^!" - 1)+. 

By Theorem 10.4 in [3], for some p > and r < 1/4 we have the estimate, 

\{v A > W) r\B 2r \< Cr n (9 + Cr a ) p {te)- p 

Let us choose r (depending on 8 to be chosen later) so that Cr a < 8. So that we have 

\{v A >t8}C\B r \< Cr n t- p = ct- p \B r \. 

Thus, by choosing t large, we will be able to make the measure of the set {va > tff) R B r a small 
factor of | Pi | independently of 8. Note that va > t0 is equivalent to wa <1 —tO. 

Let G = {v A < td} n B r . We know that \G\ > (1 - ct^ p )\B r \. The set G is also the set where 
Wa > 1 — td. On the other hand, since G C Pi, P < 1 in G, then P — wa < 40 in G. This allows 
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us to estimate the difference between —N(x) and w A c in G, where A c is the complement of the 
set A. 

Clearly w A + w A c = P-N, then N + w A a = P - w A < tO in G. Since N(x) > \/AP(x) - e 1; 
we have that in G 

w A c(x) < -N(x) +te 

< -j(l-t9)+td + Sl 

< — — if 9 and E\ are small enough (depending on t). 

Consequently, \{w A , < n B r \ > (1 - ct- p )\B r \. 

For some small n > 0, we define w c = (iD^/crx) + ^-) + - We know M + v c > —E\ in B 2 , thus we 
can apply Theorem 15. II to v c (ktx) for some small R > and get 



« c (0) <c £l + c I ^4 dy 



< Ce x + C / 1 , , , dy + C 1 , , , dy 



using that \{v c > 0} n S K -i | < Ct^n" 7 ' 



< C £l + cr/c-t"- + C I dy 



since r < 1, we can bound the third term independently of r, 



y|> 



So we can choose k and E\ so that Cei + Cn a < A/(8A) and then t such that Cn n t p < A/ (8A). 
Therefore we got the following estimate 

But this means that Wk (0) < — ^ which is a contradiction since uijf (0) = 0. 

The contradiction came from saying that P(xq) > (1 — 6*) for some xq in B r . Thus P < (1 — 9) 
in B r . In the original scale, this means that P < (1 — in B r k+i, which finishes the inductive 

step and the proof. □ 



Using Lemma [9~2l we can finally prove Theorem ll.il 

Proof of Theorem 11.11 As it was mentioned before, the case a < 1 is already covered in [3] , so 
we prove the case a € (1, 2) only. 

Let us consider the fractional laplacian of order a 

-{-AY'Mx)=c a {2-a) f 6u(x,y)-^dy 
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where the constant c a remains bounded below and above for a G (1,2). 

Let b be a bump function as in Lemma P77T1 We rewrite the fractional laplacian as 



-{-&y' 2 u{x)= / 5u(x,y)b(y)^§-p-dy+ [ Su(x, y)(l - b{y)) C ^° ] Ay 



\y\ n+a V \v\ n+a 



= w(x) + g{x) 



By Lemma W% \w(x) — w(0)\ < P(x) + N(x) < C \\u\\ Lao \x\ a for some a > and C universal 
constants. 

Since u G then g G C°° . In particular g is C a . For a universal constant C, 

\g(x)-9(0)\<C\\u\\ Loo(Rn) \x\ a 
Adding the two estimates above we get that 

\{-Ay'*u{x) - (-A)-/ 2 U (0)| < C>|| ic „ (Rn) |a:|°. 

Therefore by a standard translation of the estimate we obtained that (— A) <T / 2 it e C a (Bi/2), 
with the estimate 

(-A) ff/2 « rafR <C\\u\\ LOO{Rn) . 

But if (-A) ct / 2 m G C Q . It is a classical result that this implies a corresponding estimate for 
u G (7 CT + a (see for example PS])- So we finish the proof. □ 

Remark 9.3. We have not used the homogeneity of / in any proof in this paper. With the same 
proof we can obtain the same regularity result for equations of the form 



inf L a u{x) = inf / (u(x + y) + u(x - y) - 2u{x))K a {y) dy = + b a 

aeA aeAJ Rn 

for a bounded family of real numbers b a . This is the general form of a concave uniformly elliptic 
nonlocal operator of order a. 

For the estimates, we would have to include the values of b a in the right hand side: 

M|c" + "(B 1/2 ) < C(IM| L oo +SUp& Q ) 



Remark 9.4. The assumption u G L°°(R") is not sharp. It could easily be replaced in all 
estimates by u G i 1 (K™, l/(l + |y| CT )). We kept the L°° norm for simplicity of the exposition. 
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